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Abstract 

Context. Here we investigate the dynamical evolution of the reconnection process at an initially linear 3D null point that 
is stressed by a localised shear motion across the spine axis. The difference to previous investigations is that the fan plane is 
not rotationally symmetric and this allows for different behaviours depending on the alignment of the fan plane relative to the 
imposed driver direction. 

Aims. The aim is to show how the current accumulation and the associated reconnection process at the non-axisymmetric null 
depends on the relative orientation between the driver imposed stress across the spine axis of the null and the main eigenvector 
direction in the fan plane. 

Methods. The time evolution of the 3D null point is investigated solving the 3D non-ideal MHD equations numerically in a 
Cartesian box. The magnetic field is frozen to the boundaries and the boundary velocity is only non-zero where the imposed 
driving is for stressing the system is applied. 

Results.The current accumulation is found to be along the direction of the fan eigenvector associated with the smallest eigenvalue 
until the direction of the driver is almost parallel to this eigenvector. When the driving velocity is parallel to the weak eigenvector 
and has an impulsive temporal profile the null only has a weak collapse forming only a weak current layer. However, when the 
null point is stressed continuously boundary effects dominates the current accumulation. 

Conclusions. There is a clear relation between the orientation of the current concentration and the direction of the fan eigenvector 
corresponding to the small eigenvalue. This shows that the structure of the magnetic field is the most important in determining 
where current is going to accumulate when a single 3D null point is perturbed by a simple shear motion across the spine axis. 
As the angle between the driving direction and the strong eigenvector direction increases, the current that accumulates at the null 
becomes progressively weaker. 
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O X Introduction 

\ I Magnetic reconnection is one of the main processes responsible for releasing free magnetic energy in complicated stressed 
J> magnetic field. A better knowledge of this process therefore helps us to better understand the possible non-ideal dynamical 
k> evolution of such magnetic fields. One such example is the magnetic field filling the photosphere - corona region. It has a highly 
j_] complicated structure with variations on all possible length scales. It evolves dynamically, both through slow ideal changes 
Ci that maintain the topology of the magnetic field and sometimes through fast, local, non-ideal processes where the stressed 
magnetic field releases a significant fraction of its free energy. Astrophysical magnetic fields are three-dimensional (3D) and 
therefore to understand the reconnection process, we have to investigate generic 3D magnetic fields. In 3D space magnetic 



tliererore to understand trie reconnection process, we riave to investigate generic iu magnetic iields. in iu space magnetic 
reconnection can take place under two generally different circumstances (^ Schindler et al. 19881) . either in regions where the 



magnetic field vector is not vanishing, or at 3D null points and their associated separators. The first type of reconnection occurs 
in magnetic field regions associated with rapid, but continuous, changes in field line mapping. These regions typically allow 
for small field perturbations to generate strong current concentrations far away along the perturbed field lines, and are typically 
identified as being as sociated with Quasi-Separatrix Layers (QSLs) (Priest & Demoulin 1995). This has been disputed in more 
recent investigations ( Wilmot- Smith et al. 2010l: Pontin et al. 201 1 ) where they show that QSLs and current accumulations do 



not necessarily represent the same magnetic field line region. Reconnection may also take place in magnetic fields with nulls, 
but not only at the null points. When no null reconnection takes place in these case s it i s found to do so along magnetic 
separator lines connecting two null points ( Galsgaard & Nordlundlll996l; |Parnelletal.H2O10l) . Fi nally, reconnection ca n also 



proceed directly at magnetic null points through a number of diff'erent typical evolution processes ('Priest & Pontin"2009'). Such 
nulls have recently be shown to be present in abundance in the solar corona (Regnier et al. 2008; Longcope & Parnell 2009t). 
Aspects of rnagnetic reconnection at a siiigle rotationally symmetric 3D null point have been investigated by a number of a uthors 
jPontin et al.ll2007bHPontin & Galsgaardl l2007l:lGarsgaard & Pontinll201 lHoSsgaard et alJ[2003l:lMcLaughUn et al.ll2oTol) . Here 



the spine axis is typically perturbed by external stress and the null reacts by generating a current concentration in its vicinity 
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Figure 1 : The structure of the initial potential 3D null point for 0-0. The x-dkection represents the spine axis, the y, z-axes the 
eigenvector direction in the fan plane. Here the z-direction represents the numerically smaller eigenvalue. 



leading to a typical sp i ne-fan reconnection scenario (" Priest & Pon"tinll2009t iPontin et al.ll2007al) . To extend these investigations, 
Al-Hachami & Pontin (I2OIOI) looked at how current accumulates at a null that is not rotationally symmetric. They found the 
degree of asymmetry to have strong implications for the geometry of the current accumulation. As the asymmetry of the null 
increases the current sheet extends steadily further along the eigenvector of the numerically smallest ('weak') fan eigenvalue 
(hereafter referred to as the 'WE' direction). In their investigation, the stress of the spine axis was directed along the dominating 
fan eigenvector direction, and therefore gives no indication to the collapse of the null and the formation of the current sheet as 
the orientation betwe en the dominatin g fan eigenvector direction (SE) and the imposed stress of the spine axis changes. In a 
different experiment, iGalsgaard et al.l (11997.) investigated the behaviour of a double null pair, showing how the presence of a 
separator connecting the two null points has a clear local focusing effect on the current independent of the relative orientation 
of the imposed driver and the separator. At the same time it was found that on the opposite side to the separator, the current 
accumulation changed to become perpendicular to the imposed driving direction. Despite the difference in the topological 
structures between a double null point pair and a single null, this strongly suggests that the asymmetric structure of a null will 
be able to significantly influence the alignment of the current accumulation. This is the topic investigated in this paper A simple 
potential 3D null point is chosen, where a fixed ratio between the three eigenvalues (all different) is chosen. The influence of 
the relative orientation of the fan eigenvectors and the imposed driving on the dynamical response of the null point structure is 
investigated. 

The layout of the paper is as follows. In Section |2] we discuss the structure of the magnetic field and the imposed driving 
profiles in relation to the numerical domain. Sections |3] and |4] describe the results of the different experiments, combining both 
impulsive driving and a long time consistent stressing of the system. Section |5] discusses the implications of the findings, and 
finally Section|6]summaries the conclusions of the investigation. 



2 Numerical setup 

The magnetic field configuration adopted here is the linear rotationally asymmetric potential magnetic field defined by 

B = Bo[-2x, (1) 
zcos(6i) sin(6l)(2 - 2k) + y((2 - k) cos(ef +Ksm(0f), 
z(k cos(0f + (2-k) sm(0f) + y(2 - 2k) cos(0) sm(0)], 

where Bq is a scaling parameter, here of order unity, a- is a number between and unity determining the asymmetry of the null 
point - here we take k = 0.5 - and is the angle between the y-axis and the major fan eigenvector. (The eigenvalues for the 
eigenvectors defining the fan plane of the null are {{2-k), k).). The null point is located at {x,y,z) - (0,0,0) with the spine axis 
along the .jc-axis and the fan plane represented by the y - z plane, i.e. x = 0. The structure of the null point for 6* = is shown in 
Fig-Ill To investigate the dynamical evolution of the 3D null, the domain of interest is limited to (x,y,z) - ±(0.5, 1.5, 1.5), and 
has been discretised using a uniform grid with 200x300^ grid-points. 




Figure 2: A surface plot of the imposed normalised driving velocity profile. 



The plasma is initially assumed to be at rest and have a constant density and thermal energy. In the present experiment the 
MHD equations are non-dimensionalised and the density is set to unity, while the thermal energy is set to 0.025. 

The linear nature of the magnetic null point forces a limitation on the numerical domain size. As a Cartesian domain does 
not represent very well the asymmetry of the magnetic null structure as it is rotated around the x-axis, and because the imposed 
boundary stressing also breaks any symmetry in the 3D null point, it is a non-trivial problem to describe the boundary conditions 
as the perturbation reaches the boundaries. In this approach, the flow perpendicular to all of the boundaries is imposed to be zero, 
and parallel velocities at the boundaries are restrict to be non-zero only in the two regions where dedicated stressing velocity 
flows are imposed. 

The driving of the system is limited to two areas on the jc-boundaries, by advecting the magnetic field in the y-direction, with 
an additional spatial dependence on the z-direction. The imposed driver has the following mathematical expression, 

Vy(y,z,t) = +0.5VAt)itanh({y-yo)/yh)-tanhi(y + yo)/yh)) (2) 
X (tanh((z - zo)/zh) - tanh((z + zo)/zh)), 

where Vd{t) is a time dependent profile related to the specific experiment discussed below, and {y,z) = (iyo, ±Zq) represent the 
locations (at the centre of the hyperbolic tangent function) where the driving velocity has half of its peak value. The variables 
yi,, Zh are the half widths of the hyperbolic tangent functions. This gives a near constant driving velocity between {y,z) - 
i+iyo — yi,), ±(zo - Zh)) while the velocity decreases outside this region, asymptotically approaching zero. For these experiments 
we use yo - 0.6, yi, = 0.2, zo = 0.3 and Zh - 0.2. This profile is imposed on the two .i-boundaries with opposite directions. Two 
types of driving are used. In the first the driving is switched on at the initiation of the experiment and reaches the peak velocity 
exponentially with a given short time scale. The driving is then maintained at this amplitude for the full extent of the experiment. 
For the second series of experiments, the same spatial driving profile is used, but this time is only imposed for a limited time to 
give instead an impulsive perturbation of the system. For both cases the structure of the driving profile relative to the domain 
extent in the y - z directions is seen in Fig. [2] 



The experiments are conducted using the 3D non-ideal MHD code by iNordlund & GalsgaardI (Il997h . This is a high order 
finite diff'erence code using staggered grids to maintain conservation of the physical quantities. The interpolation operators are 
fifth order in space while the derivative operators are sixth order The solution is advanced in time using a third order explicit 
predictor-corrector method. A constant 77 model for the resistivity is assumed, with a non-dimensional value of 10"^. 



3 Continued driving 

This section discusses the case of the continuous driving. For this case the driving amplitude has been modulated in the following 
way to; 

Vrf = vo tanh(f /f„s,), (3) 

to smoothly ramp up the driving velocity to a constant amplitude. In these experiments a ramp up time of frijc =0.1 has been 
used. A number of experiments are conducted in which the angle between the direction of the driver and the orientation of the 



Table 1 : Experiments investigating the influence of the current accumulation and reconnection as the angle between the driver 
and the null point's symmetry axis is changed. 



Name ^ 


TTT^ 


Peak / 


Peak V " 




A 





36.6 


0.70 


0.019 


B 


30 


34.5 


0.62 


0.017 


C 


60 


28.7 


0.43 




D 


70 


25.4 


0.40 




E 


80 


20.5 


0.40 




F 


90 


20.9 


0.43 




lA 





6.4 


0.078 


0.000855 


IB 


15 


6.3 


0.076 


0.000755 


IC 


30 


5.9 


0.071 


0.000695 


ID 


45 


5.2 


0.061 


0.000615 


IE 


60 


4.3 


0.047 


0.000560 


IF 


75 


3.2 


0.039 


0.000465 


IG 


90 


2.7 


0.038 


0.000390 



^ First five represent continued stressing (A-F) while the last seven (lA-IG) represents the impulsive driving. 
Angle between the weak fan eigenvalue and the z-axis. 
Peak current value. 
Peak outflow velocity. 
Peak integrated parallel electric field. 



main fan eigenvector is changed between to 90 degrees, experiments A to F. To quantify the main changes to the dynamical 
evolution of the magnetic field close to the null point a number of measurable numbers are chosen. These are listed in Table [T] 
and a discussion of each of them is provided below. To simplify our discussion we define here two orthogonal directions in the 
yz-plane (the fan plane). The first is the direction of the eigenvector associated with the strongest eigenvalue, denoted the 'SE' 
direction, while we denote the direction of the eigenvector associated with the weak eigenvalue by 'WE'. 

As the driver is initiated, a pulse propagates along the field lines towards the null region. Due to couplings between the 
magneto-acoustic waves the pulse does not propagate as a pure Alfven wave and it can therefore concentrate on the null point, 
eventually creating a local current concentration at the null. As the system is continuously stressed the current magnitude 
increases and eventually degenerates the null forming a localised current concentration, at least for the majority of the conducted 
experiments. 

3.1 Current accumulation 

The frames in Fig. [3] show how the spatial orientation of the current sheet changes for different values. For = the current 
sheet is clearly tilted in the jcy-plane (the plane perpendicular to the shear driving), while extending a comparable distance in 
the z-direction, such that it approximates a tilted circular disk. As increases the orientation of the current sheet changes, 
always aligning itself along the WE direction. At the same time the width (in the SE direction) of the current sheet decreases 
making it appear increasing long. When is larger then 70 degrees, implying that the WE direction becomes almost parallel 
with the imposed driver, no significant current accumulates at the null point. Instead the strongest current is located close to the 
z-boundaries. This is most easily seen in the middle frame. The lower frame shows a transition in the current accumulation in the 
yz plane. For = the current connects in a straight line between the two z-boundaries, while for increasing it starts connecting 
with secondary current structures forming an S -like structure. This structure vanishes for experiment E and F. The effect of the 
driving is therefore different for these experiments where the driver is closest to being aligned with the WE direction. We return 
to this in the discussion of the velocity flows below. 

Fig. |4] shows the time dependence of the peak current in the A-F experiments, measured in a small region around the initial 
null position. This shows that typically the peak current magnitude slowly builds, reaching its peak value around t=7.5. The 
peak current decreases with increasing angle, 0, a similar behaviour is seen for the time where the different experiments reach 
their peak current for experiments A - D. Again experiments E and F behave differently, reaching the peak value later without 
decreasing again. This is because the current gradually extends inwards from its maxima at the z-boundaries, eventually entering 
the measurement region around the null. 




Figure 3: Images of the current structure for experiments A-F in Table[T]at t=7.5. The top frame shows the xy-plane, the middle 
frame represents the xz-plane and finally the lower frame shows the yz-plane. In all frames the redundant coordinate is zero 
and the initial null point is therefore located at the centre of the images. The dynamical range in each of the frames are scaled 
according to the absolute peak current of the experiments. 

3.2 Velocity flow 

Fig.|5]shows the time evolution of the peak velocity of the plasma, measured in a small region around the initial null position. The 
velocity has a similar dependence on as the peak current discussed above, with the peak velocity decreasing with increasing 0. 
The traditional delay between peak current and peak velocity is seen for all experiments, representing the time delay required for 
accelerating the plasma out of the diffusion region. The relative increase between experiments is larger for the velocity than for 
the current, indicating that the effective combination of the Lorentz force and the gas pressure decreases with increasing angle. 
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Figure 4: The time evolution of the peak current for the six continuously driven experiments. The lines represent experiment: A: 
full line, B: short dashed line, C: dot-dashed line, D: long dashed, E: dotted, F: triple dotted dashed line. 




Figure 5: The time evolution of the peak jet velocity for the six experiments. The lines represents experiment: A: full line, B: 
short dashed line, C: dot-dashed line, D: long dashed, E: dotted, F: triple dotted dashed line. 



Looking at the graph, it is seen that the peak velocities in experiments D, E and F follow nearly the same evolution, reaching 
a peak value later then the three other experiments. This is in contrast to the peak current, where only E and F had similar 
evolutions. 

The frames in Fig.|6]show the structure of the velocity flow around the null point for four of the experiments (A, B, and C at 
t=7.5) and {F at t=10.0) together with an isosurface of the strong current. The top frame represents experiment A, where it is seen 
that the velocity flow advects the plasma and magnetic field into the current sheet along a direction normal to the surface of the 
disk-shaped sheet, before expelling it from the narrow sides of the sheet close to the y-direction, in two high velocity reconnection 
jets. The figure shows that the current sheet is only slightly extended in the z-direction and therefore has the approximate shape 
of a tilted circular disk, but that the reconnection jets are primarily in the y-direction with only a very small velocity component in 
the z-direction. In other words for this set-up the structure of the velocity is similar to the rotationa lly symmetric null point cases 
previously investigated ("Pont in & Galsgaardll2007 : Pontin et"al]|2007at Galsgaard & Pontirill201 1 ). The second frame represents 
experiment B. Here, both the current sheet and the reconnection jets are rotated relative to the imposed driver direction {y- 
direction). The velocity flow propagates more radially away from the null providing a more homogeneous flow structure. Ahead 
of the jet, the flows are seen to align with the imposed driving direction, showing how the initial slow perturbation is imposed 
on the system, but is then subsequently distorted as the current sheet aligns with the WE direction. A similar picture is seen for 
experiment C (third frame), where it is also clear that there are differences in the alignment of the current sheet and the outflow 
velocity. The current is still aligned with the WE direction, while the tension forces are strongly influenced by the orientation of 
the imposed driving, leading to an asymmetry in the overall structure. Finally the bottom frame shows the same for experiment 
F, but for the later time, t=10. This is the time where the current reaches its peak value in a limited region around the null. But, 
as it is seen in this figure, the current accumulates inwards from the z-boundaries towards the null point as time progresses, rather 




Figure 6: The four frames show the structure of the reconnection region at t=7.5 for experiments A, B and C while F is shown 
at t=10.0. The field lines represent the velocity streamlines, with the arrows indicating the local strength of the velocity. The 
central isosurface represents the current sheet, while the two adjacent isosurfaces shows the spatial location with the highest flow 
velocities - the reconnection jet. 

than being peaked at the null as in A-C. The strong current we see for this event therefore does not represent a growing current 
sheet centred at the null. The reason for this is linked to the structure of the magnetic field. As the field is sheared in the WE 
direction, the perturbations rapidly propagate out along the SB direction reaching the z-boundary where the field line motion is 
frozen due to the imposed boundary conditions. This gives a characteristic shear motion across the fan plane building up strong 
currents close to the boundaries. For the velocity, a flow pattern is soon developed where the peak value is reached in four patches 



located around the null point lying close the x - plane. These are not associated with the jets driven by reconnection at the 
null, but are effects of enhanced diffusion away from the null due to the continuously imposed shearing and the accumulation 
of strong current close to the z-boundary. Note that the collapse around the null point is inhibited in this case because there is 
no component of the shear disturbance along the SE direction. Rather it is exactly aUgned with the WE direction, and thus the 
Lorentz forces that trigger the collapse of the null are weak owing to the weak magnetic field in the shearing plane. 

3.3 Reconnection rate 

An important parameter to measure is the reconnection rate. This gives a simple way to compare the efficiency of converting 
free magnetic energy into different types of energy for the different relative orientations. In three dimensions the reconnection 
rate is determined as the peak value of the integrated parallel electric field along all field lines passing through the diffusion 
region (Schindler et al. 1988). Here this integral in general reaches its maximum value for a field line in the fan plane. To avoid 
numerical problems integrating field lines very close to the null we integrate 

where e is a small distance from the null point — here we take a value of 0.01. L is a maximum distance to integrate to avoid 
problems from the interaction with the imposed boundary conditions. 

A simple way to check that consistent results are obtained from the experiments is to change the value of L. If the current 
is localised around the null point, then the integrated parallel electric field should converge towards a constant value as L is 
increased. When significant currents are present far from the null, the integral will continue to increasing with larger L values. 
To investigate this the integral has been calculated for 5 different values of L {- 0.5, 0.75, 1.0, 1.25 and 1.4). Fig. |7] shows the 
development of for three of the experiments (A, C, F), where the different lines represent the different integration lengths 
L. We see that only for the A and B cases do we find a convergence of ^,„a_^ for increasing L. For the other cases we found a 
continuously growing value of the parallel electric field. This shows that the continuous driving experiment is strongly limited 
in deriving reconnection rates as a function of 9. This limitation is related to the linear structure of the magnetic field and the 
imposed boundary conditions. 



4 Impulsive driving 

In this section we consider the case where the boundary driving is imposed in an impulsive manner That is, the boundary 
shearing is applied for some fixed period of time before being turned back to zero. We apply a driver with the spatial profile 
given in Eq. (|2]i but with the amplitude now being localised in time, via 

Vdit) = vo ((^f - 1) < f < 2t. (5) 

We take vq - 0.01 and r = 1.5 so that the driving switches off at f - 3. All of the other parameters are the same as in the 
continuously driven simulations described above, except that in this case to resolve the current sheet we only require to use a grid 
resolution of 96 x 128^. 



The evolution of the system is qualitatively similar to that described by lPontin et al.l (l2007al) who investigated the rotationally 
symmetric case. As with the continuously driven case described above the impulsive shear driving results in a magnetic field 
disturbance that propagates into the domain, focussing around the null point. The difference here is that for all cases a current 
layer forms in this region around the null, with the current maximum located at the null point. The magnitude and dimensions of 
the current layer increase until the driving is switched off, after which the resistive diffusion gradually dissipates the current and 
the magnetic field relaxes back towards the orthogonal null point configuration of the initial field. For longer runs oscillations in 
both angle and current magnitude were observed, as the initial perturbation was bounced between the closed domain boundaries 
and thereby re-entered the area around the null point several times. 

Here we have repeated the impulsive driving of the system with the same initial conditions as described in Section |2] The 
simulation was run with - 0,15, 30, 45, 60, 75, 90. Broadly speaking the results concerning the geometry of the current layer 
are similar to those found for the continually driven system. In particular, the extension of the current layer on large scales far 
away from the null is selected by the direction of the SE. However, locally the current flows through the null orthogonal to this, 
in the WE direction. 

The only simulation runs for which the above does not adequately describe the current sheet geometry are when 0-0° and 
9 = 90°. For these simulations the shear driving is exactly along the SE (0 = 0°) or WE {6 = 90°) directions. In each case the 
current vector at the null is perpendicular to the plane of the shear, with the current layer extending along this direction. The 
magnitude of the peak current density, however, is very different in these two simulations. Specifically, when 6 = 0° and the 
driving is along the SE direction, the collapse is much more efficient due to the higher Lorentz forces, and the peak current is 
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Figure 7: The peak integrated parallel electric field for experiments A (top), C (middle) and F (bottom) for different values of the 
integration distance, L, as a function of time. The lines represents L: 0.5 is full line, 0.75 is dotted line, 1.0 is long dashed line, 
1.25 is dot-dashed line and 1.4 is dashed triple dotted line. 



much higher. Indeed, there is a continuum of values for the peak current for the different runs, with the peak current decreasing 
as is increased, as shown in Fig. [8] The same pattern in seen in the outflow jet velocity (see Table[T]i, with the peak jet velocity 
decreasing as 6 increases and the current density decreases. 

Turning now to the reconnection rate, we again observe a similar behaviour as for the peak current density. As shown in the 
left frame of Fig.|9l the reconnection rate increases steadily as the null is sheared, reaching a maximum value around the time 
that the driving ceases (f ^ 3), before then decreasing. Examining the dependence on the angle 9 of the maximum reconnection 
rate, we find the peak value of the reconnection rate over time decreases for increasing 9 (right frame of Fig.|9|. 
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Figure 8: Left: evolution of the spatial maximum of the current density with time for the different impulsive driving runs {0 — Q 
is the soHd line, 0=15 dotted 6* = 30 dashed, 6 - A5 dot-dashed, 6 - 6Q tiiple-dot-dashed, 0-15 long dashed, 6-90 stars). 
Right: spatial and temporal maximum of |J| as a function of 9. 




Figure 9: Left: evolution of the spatial maximum of the integrated parallel electric field (i.e. the reconnection rate) with time 
for the different impulsive driving runs (0 = is the solid line, 0=15 dotted = 30 dashed, 6 - A5 dot-dashed, = 60 
triple-dot-dashed, 6-15 long dashed, 9 -9Q stars). Right: spatial and temporal maximum of the reconnection rate as a function 
of 61. 

5 Discussion 

Previous investigations of rotationally symmetric null points found the current to accumulate in the direction perpendicular to the 
imposed stress of the spine axis. This is expected and is supported using simple arguments, based on Ampere's law. When the 
rotational symmetry of the initial magnetic null point is broken, and the WE is perpendicular to the direction of the spine stress, 
the accumulating current she et extends steadily further along the WE direction as this eigenvalue decreases in relative magnitude 
( Al-Hachami & Pontinll2010l) . In this paper we showed that the current continues to accumulate along the WE direction even for 
cases where the driver direction is no longer perpendicular to this plane. This can be understood as follows: It is the Lorentz 
force that is responsible for initiating the collapse of the magnetic field around the null to form the current layer (while the plasma 
pressure force acts to oppose the collapse). Thus, so long as some component of the shear disturbance acts in the SE direction, 
the magnetic field will preferentially collapse in this plane since the Lorentz force will be strongest in this plane. The collapse 
of the null in the SE plane naturally leads to a current vector perpendicular to that plane, i.e. along the WE direction. However, 
this collapse to form a current layer at the null does not always dominate the dynamics. When the angle between the SE and 
the driving direction reaches about 70° there is a change in the behaviour of the current accumulation found for the continuous 
driving case. For these situations the imposed stress propagates out along the field lines away from the null in a direction that 
is near perpendicular to the imposed stress, and instead it accumulates close to the z-boundaries. This effect is enhanced by the 
imposed boundary conditions, that freeze the field line footpoints at the side boundaries preventing the stress from propagating 
away from the null region. As a result of the two competing effects, the current layer has a 'zig-zag' shape when viewed along 
the x-direction looking down on the fan plane, see Fig. [3] This shows that having an asymmetric 3D null and imposing a shear 
perturbation, a strong current will not necessary accumulate at the null point, depending on the orientation of the shearing with 
respect to the null structure. 



Comparing the continuous and impulsive driving experiments, we see that it is only when the system is driven continuously 
that the current begins to accumulate away from the null. In the impulsive driving case the current is always focussed at the null, 
as there is no time for the information to communicate with the domain boundaries before the pulse has passed the null point 
region. One therefore does not expect to see this domain effect for this type of driving. This is also seen in the fact that the 
peak current reached in the two experiments differs with a factor of 6, with the systematic driving being able to impose a much 
stronger deformation of the magnetic field around the null. 

Magnetic reconnection is often associated with fast plasma outflows. Here we find that knowing the structure of the magnetic 
field in the vicinity of the null point provides insight as to which directions one can expect to find the reconnection jets. The 
current aligns locally along the WE direction, with the reconnection jet being ejected mainly perpendicular to this along the SE 
direction. As the angle between the SE and the driver reaches near perpendicular, the null point is stressed in a way where no 
jet is formed due to a current concentration around the null, but the reconnection speed is still very high. Instead four jet-like 
structures are located in the old fan plane and are due to the influence from the imposed boundary conditions. This is important 
in that it is not obvious that signatures of fast magnetic reconnection in relation with nuU reconnection will always show up as 
narrow single jet regions around a common diffusion region. 

The evolution of the reconnection rate is not so 'smooth' as for J (see Figs.|2l|9]l due to the difficulty in obtaining this quantity. 
This is because E\\ is typically highest in the weak field regions (see the above discussion regarding J), so that obtaining a high 
field line density in this region to accurately capture the spatial maximum is not straightforward. The results obtained show 
a clear difference between the impulsive and sustained driving. For the impulsive driving, we find that the reconnection rate 
decreases as we increase 9 mirroring the behaviour of the peak current (which is natural since = riJ\\)- In this case the effect of 
the boundaries cannot influence the dynamics at the null region before the peak of the reconnection process has occurred. On the 
other hand, for the sustained driving only the cases with small values have the current concentrated around the null such that 
actual reconnection rates for the null reconnection can be obtained. For the remaining cases we convincingly showed that the 
contribution to determination of the reconnection rate continues to the edge of the domain. To first order this is an effect of the 
imposed boundary conditions and the linear structure of the magnetic field used. It will therefore be interesting and important to 
redo this investigation with a different magnetic configuration that Umits the influence of the faraway boundary conditions on the 
local evolution at the null. By contrast to the continuously driven case, the time-limited driving gives a current pulse that always 
becomes localised around the null. The reconnection rate therefore does not depend on the integration length L along the field 
line (so long as L is taken to be sufficiently large) as was seen for the continuous driving case. 



6 Conclusion 

Recent investigations of magnetic reconnection at single 3D null points clearly show the importance of the combination between 
the magnetic field structure and the imposed driving. While the collapse of the null - and therefore the formation of the current 
layer - depends strongly on the structure of the magnetic field, the driving seems only to be the provider of stress necessary to 
foster the collapse and the subsequent reconnection process. Both the current density in the sheet and the associated reconnection 
jet velocity are strongly related to the relative orientation between the magnetic field structure and the driver, with the cases where 
the null point's SE direction align with the stress of the spine axis providing the largest amplitudes of these directly measurable 
physical variables. 

In order to model coronal magnetic fields, one can use different techniques to extrapolate the photospheric field into the 
corona, employing potential or non-linear force free extrapolations. Depending on which approach is used the local structure of 
the magnetic field around the null will change. Therefore it is difficult to give a precise description of the dynamical evolution 
of the magnetic field that would be expected when it is exposed to external (boundary driven) stresses. However, combining 
observations for the flow direction (when possible) and photospheric stressing patterns with field extrapolations, one may be 
able to predict the local behaviour of the 3D null point using the results discussed in the recent papers on 3D null reconnection 
dPontin & Galsgaard"2Q07l: IPontin etal]|2007al: kl-Hachami & Pontinll2010l: lOalsgaard & Pontinll201 i V This is what was done 



by Masson et al. (2009), but for a case with an indirect stressing of the magnetic null that led to the collapse that facilitated the 
ongoing magnetic reconnection. A similar indirect perturbation has not been investigated for generic null points, and therefore 
direct predictions of the null evolution are uncertain. 

To provide a better working hypothesis for null point behaviour more experiments need to be conducted, using magnetic 
field configurations that are not linearly dependent on space coordinates to better see how non-linear structural changes of the 
magnetic field will influence the reconnection process. It is vital to repeat the sustained driving case studied here with a different 
initial null point configuration that limits the implications of the boundaries and allows for an independent investigation of the 
reconnection rate at the null point. Furthermore, in order facilitate more direct comparison of the results with realistic situations 
we need to make scaling predictions on the behaviour of the null evolution. This is a very difficult task as it requires that we 
repeat experiments with the magnetic resistivity changing over several orders of mag nitude - a process requiring a substantially 
higher numerical resolution than the present experiments ( Galsgaard & Ponthil 201 11) . 
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